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ABSTRACT 

We study the formation of disk-dominated galaxies in a ACDM universe. Their existence is considered to 
be a challenge for the ACDM cosmology, because galaxy mergers isotropize stellar disks and trigger angular 
momentum transport in gas disks, thus fostering the formation of central stellar spheroids. Here, we postulate 
that the formation of stellar spheroids from gas-rich disks is controlled by two parameters that characterize 
galaxy mergers, the mass ratio of merging dark matter halos, and the virial velocity of the larger merging 
halo. We utilize merger histories generated from realizations of the cosmological density field to calculate the 
fraction of dark matter halos that have avoided spheroid formation, and compare the derived statistics with 
the spheroid occupation fractions in surveys of nearby galaxies. We find, for example, that the survival rate 
of disk-dominated galaxies in ACDM is just high enough to explain the observed fractional representation of 
disk-dominated galaxies in the universe if the only mergers which lead to central spheroid formation are those 
with mass ratios Mi/M\ > 0.3 and virial velocities V v ; r .i > 55 km s -1 . We discuss the physical origin of this 
criterion. 

Subject headings: cosmology: theory — galaxies: formation — galaxies: spiral 



1. INTRODUCTION 

The existence of disk-dominated galaxies, with little 
or no bulge, is frequent ly cited as a challenge to the 
ACD M cosmology (e.g., iKautsch et"ai1 120061: iKormendvl 
120071) . Apart from the question of whether or not 
the theory of galaxy formation in ACDM can suc- 
ceed in making these observe d rotationally-supported disk 
galaxies in the first place jD' Onghia &^urkertl l2004t 
lAbadi. Navarro. & Steinmetzl l2003h . the survival of such 
disks, once formed, is our focus here. The mergers that 
every galaxy-hosting dark matter halo experiences can trig- 
ger angular momentum transport in the interstellar medium 
of the merger remnant. If a substantial amount of gas is 
transported into the central kiloparsec of the remnant, the 
gas can fuel a starburst forming a central stellar system — 
a "classical" bulge (i.e., self-gravitating, baryon-dominated 
stellar system supported by random motions rather than 
rotation). By contrast, "pseudobulges" can also arise in 
some disk galaxies if they have not experienced recent ma- 
jor mergers, by the secu lar transport of angular momen - 
tum (e.g., by galactic bars. iJo gee. Scoville. & Kenneyi 2005). 
These pseudobulges (sometimes called "disky bulges") are 
suppo rted more by rotation than r andom motion, however 
(e.g.. lKormendv & Kennicuttl l2004. and references therein). 3 
About one third of all disk galax ies in the local univer se do not 
contain bulges or pseudobulges ( Kautsc h~et al.l2006h 4 and an- 
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3 Pseudobulges are to a larger degree supported by rotation than classi- 
cal bulges and can have rotational velocity-to-lD velocity dispersion ratios 
Vmax/o" > 1. Their velocity dispersions are smaller than expected from the 
Faber-Jackson relation. Pseudobulges tend to have flattened isophotes and 
surface brightness profiles close to exponential. The classical bulges and 
pseudobulges can be distinguished from a third class, the boxy or peanut- 
shaped bulges, wh ich are bars seen edge-on (Bureau & Freemaij 119991 ; 
Athanassoula 2005). 

4 The classification of galaxies by their bulge content is usually done using 
simple bulge-disk or concentration m odels, which may not be accurate in e x- 
treme disk-dominated galaxies (e.g., Boker, Stanek, & van der Marel l2003h . 



other third contain only pseudobulges. lAllen et all d2006l) car- 
ried out a Sersic spheroid and exponential disk decomposition 
on a large sample of galaxies and find that 30% of exponential 
disks h ave small bulge-to-total ratio B/T < 0.2. Barazz a"et alJ 
(2007) report that 20% of disk galaxies can be visually classi- 
fied as bulgeless. 5 

For our purposes here, we shall assume that the survival 
of disk-dominance means that no classical bulge is formed. 
Although the precise characteristics of mergers that form 
classical bulges remain unknown, disk-dominated galaxies 
must have avoided major mergers during and after forma- 
tion. Just how large the mass ratio of the merging galactic 
halos must be in order to induce bulge formation is some- 
what uncertain, however. Mergers of similar-mass galaxies 
have been shown to trigger starbursts and result in ellipti- 
cal galaxies, by gas dynamical and TV-body simulations for 
galaxy halo masses M > 10 n M Q , mergi ng at relative veloc- 
ities o f the order of th e virial velocity dMihos & Hernquisa 
Il996t ICox et al.1 120071) . In that case, a bulge forms from 
the momentum-exchange and compression of gas in merger 
shocks and the outward angular momentum transport induced 
by merger torques. Pure TV-body simulations of such merg- 
ers find that th e pre-existing stellar d isks are mixed and 
destroyed (e.g., iNaab & Burkertl 120031) . leaving an ellipti- 
cal galaxy. For minor mergers, TV-body and gas dynami- 
cal simulations with halo masses M > 1 n M w indicate that 
disks s urvive but bulges c an also grow (Mihos & Her nauisl 
1 19941: ICoxetalJ 120071: lEliche-Moral. Balcells. & Aguerril 
I2006L see also iD'Onghia et alJ 120061) . However pure TV- 
body simulations of this process find that these bulges 
are pseudobulge-like (i.e., partially supported by rotation), 
rather than classical, for mass ra tios 0.1 < M2/M1 < 0.25 
dBournaud. Jog. & Combes 2005). Bulge formation by the 
merging of lower mass (V max < 70 km s -1 ) and gas-rich galax- 
ies has not yet been simulated. 

s Pure disk galaxies a lso contain nuclear star clusters (e.g., B oker et all 
2002; Walcher et al. 2005) which could be products of secular gas transport 
iMilosavlievic 2004), but these star clusters are tiny by comparison with the 
bulges described above. 



2 



THE ABUNDANCE OF DISK-DOMINATED GALAXIES 



Semianaly tic models of galaxy for mation 
(e.g., [K auffman nTWhite. & Guiderdonil 119931: 

Baug h. Cole. & Frenkl I1996T assume they can track the 
morphological type of galaxies by converting a disk com- 
ponent into a spheroidal component in mergers with mass 
ratios p, = M2/M1 greater than some threshold. These 
models successfully reproduce the distribution of earlier 
morphological types by tuning the critical mass ratio for 
disk destruction, and adopting a critical bulge-to-total mass 
ratio that discriminates broadly between disk and elliptical 
galaxies. Recent semianalytic models employing merger 
trees extracted f rom hydrodynamic A^-body simulations yield 
a similar result (Ma iler et al.ll2006 b. However, those studies 
which focus on the survival of disks generally do not address 
the abundance of disk-dominated galaxies. Existing ab initio 
cosm ological simulations have yielded disks with bulges 
(e.g. jAbadi. Navarro. & Steinmetzll2003l) . but currently lack 
the dynamic range to explore a large enough volume to sam- 
ple the statistics of galaxy morphology while simultaneously 
resolving the formation and mergers of individual galaxies. 

The purpose of this work is to compare the predicted disk 
survival probabilities during hierarchical merging in a ACDM 
universe with the observed statistics of galaxies at the end of 
the Hubble sequence. This comparison is then used to place 
constraints on the physics of bulge-forming mergers. In § 12.11 
we discuss the role of mergers in bulge formation. In § 12.21 
we describe the effect of cosmic reionization and Jeans-mass 
filtering on bulge formation. In § 12.31 we present a calcula- 
tion of bulge formation and disk survival probabilities derived 
from galactic halo merger trees ge nera ted from realizations of 
cosmological density fields. In § 12.41 we compare observed 
disk galaxy statistics with these merger tree results and place 
constraints on the properties of bulge-forming mergers. In 
§ |3] we summarize our main conclusions. Standard cosmo- 
logical parameters consistent with the Wilkinson Microwave 
Anisotropy Probe dSpergel et al.ll2007l) are assumed through- 
out. 

2. DISK SURVIVAL IN ACDM COSMOLOGY 
2.1. Mechanisms for Bulge Formation in Mergers 

The merging of gas-rich galaxies fosters bulge formation 
directly and indirectly. Directly, the time-dependent gravita- 
tional potential of the two merging components deflects some 
of the gas into the center of the merger remnant, where it gets 
compressed in shocks and fuels a starburst. Indirectly, the 
gravitational tidal field of the merging components excites 
nonaxisymmetric perturbations inside the merging galaxies 
(bars, spirals, etc.) w hich then torq ue disk gas into the cen- 
ter of the galaxy (e.g.,[Combes 1998, and references therein). 
In the center, again, shocks are ubiquitous and play a role 
in angular momentum transport. The indirect channel should 
be important in minor mergers, especially where the smaller 
galaxy loses its gas to ram pressure stripping in the early 
stages of the merger. 

The strength of direct merger torques is a function of the 
mass ratio of the host dark matter halos of the merging galax- 
ies, /.i = M2/M1. The strength of nonaxisymmetric distor- 
tions in minor mergers with \i <C 1 depends nontrivially on 
the resonance st ructure of the merging halos and the orbit of 
the merger (e.g., Goldreich & Tremaine 1980). These cannot 
be modeled in a general case; therefore, we here consider only 
the gross properties of the mergers, averaged over the merger 
orbital parameters and over the properties of merging halos 
with given masses at a specific redshift. 



Efficient transport of angular momentum in perturbed gas 
disks requires an excitation of nonlinear waves by nonaxisym- 
metric potential distortions that can give rise to momentum- 
transporting shocks. 6 The strength of merger shocks is char- 
acterized by the Mach number M S h, which is the ratio of the 
shock velocity V s h to the sound speed of the warm neutral gas 
c s ~ lOkms -1 . Merger-driven strong shocks are radiative, 
because the post-shock cooling time is much shorter than the 
dynamical and sound crossing times of the H I disk. In this 
limit the shocks are isothermal and the shock compression is 
-Ml. 

In major mergers, we expect V s h ~ V ga i, where V ga i is the 
relative velocity of the two galaxies, while in minor mergers 
V S h ^ Vg a i, although the forcing of the gas can be strong where 
the gravitational torque is amplified locally by a resonance. 
In view of these considerations, we postulate that, besides 
the mass ratio fi, the efficiency of gas transport in mergers 
is controlled by a second parameter, the merger Mach number 
Ai mei = Kir,i/c s ~ /3 _1 Vgai/c s , where V v j,-,i is the virial veloc- 
ity of the larger merging halo, and j3 ~ 1 is a dimensionless 
ratio of the orbital velocity V ga i of the two galaxies (i.e., merg- 
ing dark matter halo centers) at the small radii that are relevant 
to bulge formation (r ~ few kpc) to the virial velocity Vvir,i. 7 

The relative velocity of gas disks could be smaller than 
the virial velocity, i.e., j3 < 1, if the merger starts affect- 
ing the gas only when the distance between the halo cen- 
ters has become much smaller than the scale radii r s in the 
iNavarro. Frenk. & White! d 19971 NFW) profiles of the two ha- 
los. The circular velocity of the NFW dark matter density 
profile p oc (r/r s )"'(l + r/r s )~ 2 at small radii r <C r s equals 
Vci rc ~ 1-3 Kir(>"/ r s) 1 ' /2 I where the factor 1.3 is for halo con- 
centrations c = r2ooA.s = 10. The factor depends only weakly 
on c; it is only slightly smaller (« 1.1) for c = 5. The circular 
velocity V c j rc reaches its maximum at r«2.2 r s . The scale 
radius equals 

r^lSkpcc- 1 (j^j^j ' [fUl+z) 3 + 0\r I/3 /*- 2/3 . 

' (1) 
With this we find that 

/c\i/2/ r \ 1/2 ( M \" 1/6 
(0.42+0.08.) (-) (_) , 

(2) 

implying that at low redshifts, the relative orbital velocity of 
the galaxy centers could be half the virial velocity of the larger 
merging halo when the separation is r ~ 1 kpc. 

We emphasize, however, that in minor mergers, one must 
not attempt to identify M.^ directly with (3A4 mer , because 
while A^mer 3> 1 in all cases of interest, the shock Mach num- 
ber, which reflects the velocity of secondary gas flows in a per- 
turbed merging system, can well be close to or below unity, 
implying a regime in which angular momentum transport is 
inefficient. The true condition for bulge formation triggering 
could be /i(/3-M me i-) Q > fait, where / cr i t is a threshold and a 
is a power, e.g., a = 1 (the "linear hypothesis") and a = 2 (the 
"compression ratio hypothesis"). 

6 For angular momentum transport by spiral shocks, see, e.g. 
IRozvczka & Spruitl fl993l) . [Savoniie et all fi994l) . and lGoodman & Rafiko\j 
1200 II) . 

7 Here and throughout the letter, V Ya = (GM^oo/^oo) 1 ' 2 . where M200 is 
the mass inside a sphere with radius r2oo centered on the halo within which 
the mean density equals 200 times the critical density of the universe. 
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An important distinct possibility is that in which bulge 
formation is dri ven cumulatively, ra ther than induced in a 
single merger (Bournaud et al.l 120071) . The central density 
of gas in the disk of a late-type disk galaxy could increase 
gradually due to slow, continuous radial gas inflow. Evi- 
dence for such inflow can be found in the "central light ex- 
cess" (above the exponential law) in p ure disk galaxies (e.g., 
iBoker, Stanek, & van der Marell 120031) . The inflow could 
be excited by perturbations associated with minor mergers. 
Through their differential gravitational perturbations, many 
consecutive minor mergers can induce a slow, secular drift 
in the angular momentum distribution of the disk fluid, which 
could lead to central accumulation without leaving any char- 
acteristic signatures of merger-driven evolution (dynamically 
hot stellar components, etc.). The inflow could also be driven 
by nongravitational proce sses, such as the m agnetorotational 
instability in the gas disk (Milosavlievic 2004, see also, e.g., 
Piontek & Ostriker 2004). The resulting increase of the cen- 
tral surface density brings the galaxy closer to the thresh- 
old for gravitational instability. The bulge or pseudobulge 
formation-triggering merger then must only nudge the galaxy 
over the threshold for, e.g., nuclear bar formation, where the 
galaxy is already marginally unstable. 

For systems which are gas-poor, collisionless mergers of 
stars can result in an elliptical galaxy or a c lassical bulge, if 
and only if the n is sufficiently large ( e.g., Naab & Burkert 
120031: iBournaud. Jog. & Combes! 12005). This bulge forma- 
tion can be characterized by merger ratio fi, but it is indepen- 
dent of y v ir i because the gravitational dynamics without gas 
is scale free. However, there will still be a dependence on 
Vvirj through the requirement that the galaxies prior to their 
merger were able to form long-lived (i.e. low-mass) stars (see 
below). 

2.2. The Critical Virial Velocity for Bulge Formation after 
Reionization and Jeans-Mass Filtering 

In order for the merger of two halos to have produced a 
bulge, the halos must have contained a substantial amount of 
gas, or else stars already formed from collapsed gas. The 
gas content of small-mass halos, however, was affected by 
the reheating of the intergalactic medium (IGM) out of which 
the gas inside those halos collapsed, by cosmi c reionization, 
a phe nomenon known as Jeans-mass filtering (Shapi rcTet al.l 
1994). The gas pressure of the reheated intergalactic medium 
competes with gravitational instability, in that case, to sup- 
press structure formation in those baryons which would oth- 
erwise have formed galaxies with virial velocity below some 
threshold. The Jeans length in the IGM for a gas photoheated 
to ~ 10 4 K corresponds to a halo mass after collapse and viri- 
alization for which the circular velocity is 

Vdrc = 55(r IGM /10 4 K) 1 / 2 kms- 1 (3) 

dlliev et al.ll2007l) . The actual threshold virial velocity is un- 
certain, because one must account for the time-dependent 
growth of fluctuations in an evolving background and because 
the formation of dark matter halos affects the baryons in a 
nonlinear way. Estimates of the value of the velocity threshold 
which results range from about 30 to 80 km s" 1 (|Efstathioul 
1992t iThoul & W einberg 119961: iNavarro & Steinmetzl 119971 
Kitavama & Ikeuchi 2000). Whatever the precise value 
should be, this would impose a lower limit to the critical virial 
velocity of merging halos capable of producing a bulge, as de- 
scribed above, i.e., Vcrit.mm ~ 30- 80 km s _I . 



Since the virial velocity threshold which results from Jeans- 
mass filtering depends primarily on the temperature, other 
sources of IGM heating could have a similar effect. The su- 
pernova explosions associated with massive star formation, 
for example, could also heat the intergalactic gas. Such feed- 
back could also have unbound the interstellar gas from the 
galaxies which formed these stars, if the galaxy virial veloci- 
ties were small enough. 

2.3. Merger Histories of Low-Mass Galaxies 

To explore the sensitivity of the fractional abundance of 
disk-dominated galaxies produced during structure formation 
to the critical values of /i and V v i r .i, and thus to place con- 
straints on the values of these two parameters that are com- 
patible with the observed statistics, we generate merger his- 
tories of low-mass, disk-galaxy hosting halos and study disk 
survival in this population of halos. We utilize those merger 
histories to calculate the abundance of disk-dominated galax- 
ies as a function of ^ C1 -i t and V v ii-,crit- We compare the result- 
ing abundances with the incidence of late-type galaxies in the 
Tully Galaxy Catalog (§ I2.41 >. 

The merger histories are generated from the nonlinear evo- 
lution of the initial, linearly perturbed cosmological den- 
sity field using the p ublicly-available Lag rangian perturbation 
code PINOCCHIO dMonaco et alJ |2Q02b. The code generates 
a Gaussian-random field on a cubic mesh, distributes particles 
on the mesh, and determines the collapse time of each parti- 
cle using an ellipsoidal collapse criterion. The "collapsed" 
particles are moved by Lagrangian perturbation theory and 
related to virialized objects, which are the dark matter halos, 
by a linking criterion. We employed 512 3 particles with cos- 
mological parameters f2 m = 0.24, ft\ = 0.76, a% = 0.74, and 
h = 0.73, in a cubical box with 50 comoving Mpc on a side. 
The mass of an individual particle was m part = 3.3 x 10 7 M Q , 
and halos with more than 10 particles were selected for in- 
clusion in the merger tree. For a given redshift, PINOCCHIO 
provides a list of all of the halos with mass M > 10 m part which 
formed inside the comoving box at this or any earlier redshift, 
and a complete list of their merger events. Each merger event 
is characterized by the merger redshift and the masses of the 
halos participating in the merger. 

We compute the fraction of halos containing disk- 
dominated galaxies as a function of the threshold for spheroid 
formation that is parametrized by the critical mass ratio /i 
and critical virial velocity V v i r i of the larger halo at the time 
of each merger. Specifically, we assume that a merger with 
A* > Merit will create a central stellar spheroid if the halo has 
a virial velocity VW,i > Kir.mt- We follow the most massive 
progenitor branch of the merger history of each halo and iden- 
tify the resulting z = halo as containing a disk-dominated 
galaxy if no spheroid has yet formed in the halo based on the 
defined criterion. 

There are rare cases in which a progenitor mass is so small 
at high redshift that bulge-forming mergers are not resolved 
by our numerical results. We have estimated the number of 
such cases and confirmed that it is negligible. For 1/10 < 
Merit < 1/3, the fraction of current halos in the mass range 
we will describe below, in which bulge-forming mergers oc- 
cur with a halo containing less than 50 particles is less than 
3%. For /x cr j t = 1 /2, the fraction increases to 8%, but the total 
mass of the merger remnant is small for those mergers with 
/i > 1/2. If the smaller halo contains 50 particles, the merger 
remnant with fj, > 1/2 has at most 150 particles, which is ten 
times smaller than M m ; n , the minimum mass of interest for our 
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comparisons with present-day galaxies. 

We consider halos with present masses in the range M m j n < 
M < M max , where M min = 5 x 1O 1O M (corresponding to 
Kir.min ~ 60 km s -1 ) and M max = 10 12 M Q (corresponding to 
V v ir,max ~ 160 km s~ ); the resulting galaxy statistics are com- 
pared with the observed galaxy statistics in the same approx- 
imate mass range. [Since the halo mass and the maximum 
circular velocity of the galaxy disk are not known for most of 
the galaxies in each observed sample, we use the Tully-Fisher 
relation to estimate halo masses for the observed galaxies.] 
The present total number density of halos in the above range 
is 0.021 Mpc~ 3 . Disk-dominated galaxy abundances thus cal- 
culated will not be strongly dependent on the specific choice 
of M max because halos with masses M < M max dominate the 
number density of halos in the universe t oday . However, the 
abundances will be sensitive to M m \ n . In 32.41 we explore the 
sensitivity to the choice of M m \ n . An approach more accurate 
than the one employed here would dispense with M m { n and 
would consider halos of all masses and then match the frac- 
tional disk and irregular galaxy abundances as a function of 
halo mass. 

We ignore the finite duration of the merger, which is the 
time elapsed between the halo contact and the final bulge for- 
mation. Indeed, the delay accounting for a finite merger dura- 
tion will affect only the low-redshift, disk-destroying mergers 
which are in the minority (z < 1, Fig. 0J. The dynamical fric- 
tion time scale w as recently calibrated in A^-body simulations 
(Jian g et al.ll2007h . 
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where e is the circularity parameter. Setting e = 0.5 and Vdrc = 
V_oo, this simplifies to 
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where H(z)~ l is the Hubble time at halo merger. If a pair of 
halos with mass ratio fj, > 0.1 merges before z = 1, the galax- 
ies in the halos merge by z = 0. In the halo mass range corre- 
sponding to disk-dominated galaxy hosts, the present fraction 
of mergers in progress is only < 5% for mass ratios /i > 0.2 
and< 10% for ^> 0.1. 

2.4. Results 

Figure Q] shows the fraction of galaxies without classi- 
cal bulge as a function of the bulge formation criterion 
( VWcri ti Mcrit)» which is the result of the model described in 
§ 12.31 In order to compare these theoretical contours with 
the observed abundance of disk-dominated galaxies, we select 
2281 galaxies in the nearby universe from the Tully Galaxy 
Catalog 9, that have blue magnitudes in the range -20 < Mb < 
-17 and are located at distances D < 20 h Mpc at which 
the catalog is reasonably complete. This luminosity range is 
chosen to render the number density of galaxies equal to the 
density 0.021 Mpc -3 of halos that we synthesize (see §[23}. 
The corresponding circular v elocity range calculated from the 
Tully-Fisher relation (e.g., iKannappan. Fabricant. & Franxl 
2002) is 60 km s -1 <V C < 160 km s _I . This range is consis- 
tent with the range of virial velocities in our theoretical halo 
sample, which is a self-consistency check of our association 

8 http://haydenplanetarium.org/universe/duguide/exgg_tully.php 
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FIG. 1 . — The fraction of disk-dominated galaxies (galaxies without classi- 
cal bulges) that results from bulge formation criteria (Vvir.cribA'crit) character- 
ized by the critical merger mass ratio, fi C n t , and the critical virial velocity of 
the larger halo at merger, V v ; r cr ; t . 
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FIG. 2. — Same as in Fig.fT] but the contours are labeled by the cumulative 
fraction of morphological types in the Tully sample of galaxies. Assum- 
ing that the type Sc and later do not contain classical bulges, the formation 
criteria along the Sc contour are compatible with the observed fraction of 
disk-dominated galaxies in the sample. 



of galaxies with dark matter halos in our numerical halo cata- 
log for ACDM. 

The assumption that the Tully-Fisher relation can be used to 
relate the luminosities of galaxies in the Tully catalog to their 
halo virial velocities (and, hence, to their masses, M200) works 
best for the spiral galaxies but is less certain for the elliptical 
and SO galaxies. The Tully-Fisher relation for SO galaxies is 
shifted to lower luminosity by about Mb ~ +1. 5 for each V m 



and t he scatter is larger compared to spirals dBedregal et alj 
2006). This estimate is uncertain because the Tully-Fisher re- 
lation or virial mass-to-light ratio is unknown for SOs at small 
mass near M m \ a . The virial mass-to-light ratio for elliptical 
galaxies could be a factor of 10 larger relative to spirals in the 
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fi-band dHoekstra et alj|2005t iGuzik & SeliakH2002h . If we 
shift the luminosity range of the subsample that corresponds 
to M m i n < M < M max by +5 mag for ellipticals and +1 .5 mag 
for SOs, respectively, then all the ellipticals will be removed 
(i.e., M > M max ) and the number of SOs will increase by 30%, 
but the sum of E and SO will only decrease from 25% to 20%, 
and the Sc contour move from 34% to 36%. Hence, uncertain- 
ties regarding the virial velocities of the ellipticals and SOs in 
the sample do not significantly affect our determination of the 
disk-dominated portion. 

Figure [2] is the same as Figure [T] it shows the fraction 
of disk-dominated galaxies for various bulge formation cri- 
teria. Each pair of morphological type and fraction printed on 
the contour indicates that that morphological type and later 
types occupy the corresponding fraction in the subsample of 
the Tully catalog. If we choose to assume that a particular 
morphological type (e.g., Sd, Sc, or Sb) and all later morpho- 
logical types are disk-dominated galaxies, while the earlier 
morphological types are galaxies with classical bulges, then 
the fraction of disk-dominated galaxies in the subsample is 
explained by the parameters (Vva tC nu Merit) along the contour 
labeled by the chosen transitional morphological type. The 
classification by morphological type does not precisely sepa- 
rate galaxies that have classical bulges from those do not; nev- 
ertheless we assume that Sc and later type galaxies are either 
bulgeless or have pseudobulges, while Sbc and earlier types 
contain classical bulges. Assuming this correspondence, the 
contour "Sc" in the figure shows the parameter space locus 
yielding the abundance of galaxies without classical bulges. 

While the criterion for the survival in mergers of a given 
morphological type must lie on the appropriate contour, from 
the statistics alone it cannot be determined which specific 
value of (/^crit, Kir,crit) along the contour is the true physi- 
cal criterion for bulge formation. For example, the observed 
abundance of disk-dominated galaxies is consistent with the 
hypothesis that mergers with \i > 0.3 in halos with V V a t i > 
55 km s -1 create classical bulges, while mergers that do not 
satisfy these criteria do not. This hypothesis is not unique; a 
somewhat larger /i or j t and somewhat smaller Vvi r ,crit, and vice 
versa, would be equally plausible on the basis of the statistics 
alone. Were the critical velocity for bulge formation above 
~ 65 km s _I , however, then the relative abundance of disk- 
dominated galaxies would be greater than observed for all 
mass ratios fi. In that case, there would be too few galax- 
ies with bulges relative to their observed abundance. A sim- 
ilar upper limit to V v ii-.crit results from the fact that classical 
bulges are unlikely to form from mergers that are too minor. 
If Kir,crit ~ 65, mergers only produce enough disk-dominated 
galaxies if /i cr ; t <C 1, which may be implausibly small. On 
the other hand, the critical virial velocity cannot be much less 
than the mi nim um value imposed by Jeans-mass filtering dis- 
cussed in {32.21 This means that V v ir,crit cannot be much less 
than ~ 60 km s -1 , either. To identify the true, unique crite- 
rion for bulge formation, one must resort to physical insight 
to exclude implausible, extreme criteria that are still allowed 
by the statistics. 

If we assume that /i cr it lies between 0.4 and 0.2, then the 
critical merger Mach number .Mmer,crit = Kir,crit/(10 km s _I ) 
lies between 5 and 6. When the two halo centers have ap- 
proached to within kiloparsecs of each other, the true Mach 
number of the gravitational perturbation will be reduced by 
the value of (3 ~ i. Therefore, the question of the phys- 
ical plausibility of the criterion can be rephrased: In view 
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FIG. 3. — The redshift distribution of first mergers that create classical 
bulges. We assumed /i cr ; t = 0.3 and V v ir,crit = 55 km s -1 as the critical param- 
eters for parameters of bulge formation, which is consistent with the abun- 
dance of disk-dominated galaxies in the local universe (see Fig. |2). 



of the gravitational hydrodynamics of the gas disks in merg- 
ing galaxies, is it physically plausible that a merger with 
[i ~ 0.2-0.4 and Ai mer ~ 5-6 does not trigger bulge for- 
mation, but a merger with a larger /iA4* er does lead to central 
gas inflow and bulge formation? For the ad hoc choice a = 1 
and assuming /3 ~ i, this yields a criterion /«crit(/3-A4vii crft) Q = 
/ crit ~ 0.6-1 . 

The lKautsch et al.l d2006l) sample of edge-on disk galaxies 
gives a similar disk-dominated galaxy fraction. If we assume 
the same additional fraction of early-type galaxies as in the 
Tully subsample, the fraction of Sd is 1 1 % and the fraction of 
Sc or later is 45%. 

To test the dependence on mass cutoff M m ; n , we vary 
this parameter from the fiducial cutoff at 5 x 10 10 M Q to a 
lower cutoff at 3.2 x 10 10 M Q . The number density of ha- 
los increases by 50% to 0.034 Mpc~ 3 . To compensate for 
the change in the number density of galaxies, we move the 
lower luminosity cutoff for selection from the Tully sample to 
Mb < -15.5. The observed fraction of disk-dominated galax- 
ies (Sc or later morphological type) remains unchanged at 
w 33%, but the contours representing the bulge formation cri- 
terion move by about 10% toward lower critical virial veloc- 
ities. This is because mergers with a fixed mass ratio [i tend 
to occur at smaller V v ; r j in smaller halos. The criterion, e.g., 
with (/Xcrit; Kirxrit) = (|;57 km s -1 ) on the Sc contour in Fig- 
ure|2]shifts only a small amount, to (|,50 km s -1 ), as the mass 
cutoff is lowered to 3.2 x 10 10 Mq. Any lower cutoffs than this 
are not appropriate, because at luminosities Mb < -15, disk 
galaxies gi ve way to irregulars as the most com mon galaxy 
type (e.g., Bing geli. Sandage. & Tammannll 1988b . 

In Figure [3] we plot the redshift distribution of the earli- 
est bulge-forming mergers of z = halos for /x cr i t = 0.3 and 
Kir.cnt = 55 km s -1 . The distribution is insensitive to the 
choices of and V v ir,criti as long as the two parameters 
remain on the same contour in Figure [2] This shows that 
bulge-forming mergers generally took place long after reion- 
ization was completed (i.e., Zm ^ 6). Moreover, the typ- 
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FIG. 4. — Fraction of halos that have not experienced a merger with mass 
ratio Mi I 'Mi or larger since redshift z. 



ical collapse epoch f or these merging halos was at z < 6 
(Shapiro & Iliev 2002), so their star formation was also post- 
reionization. Hence the lower limit to V v ii.crit imposed by 
Jeans-mass filtering is applicable. 

In Figure |4] we plot the fraction of halos that have ex- 
perienced a bulge-forming merger after redshift z. The fig- 
ure shows that 60% of the halos with masses in the range 
(0.5- 10) x 1O U M0 have not experienced mergers with \x > 
0.05 after z = 1, an d 30% have not experie nced mergers with 
/i > 0. 1 after z = 2. iToth & Ostrikerl d 19921) placed constraints 
on the mergers that could have taken place during the lifetime 
of a galactic disk by quantifying the role of mergers in the 
heating and thickening of the Milky Way's disk. The Galaxy 
could not have accreted more than 5% of its present mass 
during the past 5 Gyr, they found. Subsequent work refined 
the estimates of disk heating, resul ting in less stringent con- 
straints on the merger history (e.g., Velazquez & Whitel il999t 



iBenson etal.1l2004 . IKauffmann & White! d 19931) generated 
merger histories of Milky Way-sized halos using the Press- 
Schechter excursion set theory and found that in an open uni- 
verse with il m = 0.2 and SI a = 0, the abundance of disks is con- 
sistent with the Toth-Ostriker constraint. We considered halos 
with masses smaller than that of the Milky Way, which in the 
standard ACDM universe merge less frequently (70% have 
had no merger with /.i > 0.05, and 80% have had no merger 
with /i > 0. 1 since z = 1) than the halos i n earlier studies (30% 
for ji > 0.05 and 50% for ^ > 0.1 in IKauffmann & White! 
fT993h . 



3. CONCLUSIONS 

In order to explain the observed space density and fraction 
of disk-dominated galaxies within the ACDM cosmology, we 
propose a bulge-forming criterion such that only those halo 
mergers with mass ratio greater than /i cr i t ~ 0.3 and halo 
virial velocity (of the larger halo) above V v ii-.crit ~ 55 km s _I 
formed classical bulges, while other mergers did not. This 
criterion has some degeneracy between /i cr i t and V^rit, but 
Kir.crit cannot be larger than about 65 km s~ without under- 
producing the galaxy fraction with bulges, or much smaller 
than ~ 60 km s~ since Jeans-mass filtering after reionization 
inhibits such small-mass halos from acquiring and retaining 
baryons or forming stars. This bulge-forming criterion also 
gives a reasonable dimensionless condition, /ici-it/JA^me^d-it ~ 
1, for the impact of merger on the gas disk from the point 
of view of angular momentum transport. The validity of this 
bulge formation criterion needs to be confirmed by further an- 
alytic calculation or hydrodynamical simulations of mergers 
in the halo mass range V v ir,ciit ~ 60 km s~ . 
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